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Abstract 

We study the computational complexity of the Word Problem (WP) 
in free solvable groups Sr,d, where r > 2 is the rank and c? > 2 is the 
solvability class of the group. It is known that the Magnus embedding 
of Sr.d into matrices provides a polynomial time decision algorithm for 
WP in a fixed group Sr.d- Unfortunately, the degree of the polynomial 
grows together with d, so the uniform algorithm is not polynomial in d. 
In this paper we show that WP has time complexity 0(rn Iog2 n) in Sr,2, 
and 0{n'^rd) in Sr,d for d > 3. However, it turns out, that a seemingly 
close problem of computing the geodesic length of elements in £'^,2 is A''P- 
complete. We prove also that one can compute Fox derivatives of elements 
from Sr,d in time 0{n^rd), in particular one can use efficiently the Magnus 
embedding in computations with free solvable groups. Our approach is 
based on such classical tools as the Magnus embedding and Fox calculus, 
as well as, on a relatively new geometric ideas, in particular, we establish a 
direct link between Fox derivatives and geometric flows on Cayley graphs. 
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1 Introduction 

In this paper we study the computational complexity of several algorithmic 
problems related to the Word Problem (WP) in free solvable groups. Let Sr,d 
be a free solvable group of rank r > 2 and the solvability class d > 2. We 
present here a uniform decision algorithm that solves WP in time 0(rn logjn) 
in the free metabelian group Sr,2 (also denoted by Mr), and 0{n^rd) in the 
free solvable group Sr.d for d > 3, where n is the length of the input word. In 
particular, this algorithm is at most cubic in n and linear in r and d for all 
free solvable groups Sr,d- Notice, that in all previously known polynomial time 
decision algorithms for WP in Sr.d the degree of the polynomial grows together 
with d. In fact, we prove more, we show that one can compute Fox derivatives 
of elements from Sr,d in time 0{n^rd). This allows one to use efficiently the 
Magnus embedding in computations with free solvable groups. On the other 
hand, we describe geodesies in Sr,d and show that a seemingly close problem of 
finding the geodesic length of a given clement from 5^.2 is surprisingly hard - 
it is NP-complete. Our approach is based on such classical tools as the Magnus 
embedding and Fox calculus, as well as, on a relatively new (in group theory) 
geometric ideas from |12| and |53| . In particular, we establish a direct link 
between Fox derivatives and geometric flows on Cayley graphs. 

The study of algorithmic problems in free solvable groups can be traced to 
the the work [35] of Magnus, who in 1939 introduced an embedding (now called 
the Magnus embedding) of an arbitrary group of the type F/N' into a matrix 
group of a particular type with coefficients in the group ring of F/N (see section 
l2.2l below). Since WP in free abclian groups is dccidablc in polynomial time, by 
induction, this embedding immediately gives a polynomial time decision algo- 
rithm for a fixed free solvable group Sr.d- However the degree of the polynomial 
here grows together with d. 

In 1950's R. Fox introduced his free differential calculus and made the Mag- 
nus embedding much more transparent [2Ql EH |22l [23] (see also Section 12. 3p . 
Namely, besides other things, he showed that an element w € F belongs to 
N' = [N, N] if and only if all partial derivatives of w is equal to in the integer 
group ring of F/N. This reduces WP in F/N' directly to the word problem in 
F/N. In particular, it solves, by induction, WP in 5*^,^. Again, the decision 
algorithm is polynomial in a fixed group Sr.d, but the degree of the polynomial 
grows with d - not a surprise since the partial derivatives of w describe precisely 
the image of w under the Magnus embedding. 



Word Problem in free solvable groups 



3 



A few years later P. Hall proved the finiteness of all subdirect indecompos- 
able finitely generated abelian-by-nilpotent groups. This implies that all finitely 
generated abelian-by-nilpotcnt, in particular, metabelian, groups are residually 
finite. About the same time Gruenberg extended this result to arbitrary free 
solvable groups [27j. Now one can solve WP in Sr,d in the following way. Given 
w G Sr,d, as a word in the fixed set of generators, one can start two processes in 
parallel. The first one enumerates effectively all consequences of the defining re- 
lations of Sr,d in Fr (which is possible since the group is recursively presented) 
until the word w occurs, and the second one enumerates all homomorphisms 
from Sr,d into all finite symmetric groups Sn (checking if a given r-tuple of 
elements in Sn generates a solvable group of class d) until it finds one where 
the image of w is non-trivial. However, computer experiments show that the 
algorithm described above is extremely inefRcient (though its complexity is un- 
known) . 

Another shot at WP in metabelian groups comes from their linear represen- 
tations. V. Remeslennikov proved in [5S] that a finitely generated metabelian 
group (under some restrictions) is embeddable into GL(n, R) for a suitable n 
and a suitable ring R ~ Ki x . . . x A'„ which is a finite direct product of fields 
Ki. In [54], see also [55], B. Wchrfritz generalized this result to arbitrary finitely 
generated metabelian groups G. It follows that G is embeddable into a finite 
direct product of linear groups. Since WP in linear groups is polynomial time 
decidable this implies that WP in G is polynomial time decidable. Notice, that 
it is unclear if there is a uniform polynomial time decision algorithm for WP in 
arbitrary finitely generated metabelian groups. 

In comparison, observe, that there are finitely presented solvable groups of 
class 3 with undecidable WP. In [34] O. Kharlampovich constructed the first 
example of such a group by simulating a universal Minski machine in WP of the 
group. There are several results which clarify the boundary between decidability 
and undecidability of the word problems in solvable groups, we refer to a survey 
[35] for details. 

Our approach to WP in free solvable groups is based on the Fox Theorem 
mentioned above. Using binary tree search techniques and associative arrays we 
were able to compute Fox's derivatives of elements w of a free solvable group 
Sd,r in time 0{n^d), where n = \w\. Significance of this result goes beyond WP 
for these groups - it gives a fast algorithm to compute images of elements under 
the Magnus embedding. This opens up an opportunity to solve effectively other 
algorithmic problems in groups Sr^d using the classical techniques developed for 
wreath products of groups. 

In the second half of the paper. Section |4l we study algorithmic problems 
on geodesies in free metabelian groups. Let G be a group with a finite set of 
generators X = {xi, . . . ,Xr} and ^ : F{X) ^ G the canonical cpimorphism. 
For a word w in the alphabet X^^ by \w\ we denote the length of w. The 
geodesic length lx{g) of an element g ^ G relative to X is defined by 

lx{g) = min{|u)| | w £ F{X),w'^ = g}. 

We write, sometimes, lx{w) instead of Ixiw'^). A word w € F{X) is called 
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geodesic in G relative to X, if |w| = lx{w). We are interested here in the 
following two algorithmic search problems in G. 

The Geodesic Problem (GP): Given a word w £ F{X) find a word u £ 
F{X) which is geodesic in G such that = u^. 

The Geodesic Length Problem (GLP): Given a word w G F{X) find 
lx{w). 

Though GLP seems easier than GP, in practice, to solve GLP one usually 
solves GP first, and only then computes the geodesic length. It is an interesting 
question if there exists a group G and a finite set X of generators for G relative 
to which GP is strictly harder than GLP. 

As customary in complexity theory one can modify the search problem GLP 
to get the corresponding bounded decision problem (that requires only answers 
"yes" or "no"): 

The Bounded Geodesic Length Problem (BGLP): Let G be a group 
with a finite generating set X . Given a word w G F{X) and a natural number 
k determine ii lx{w) < k. 

In Section im we compare in detail the algorithmic "hardness" of the prob- 
lems WP, BGLP; GLP, and GP in a given group G. Here we would like only to 
mention that in the list of the problems above each one is Turing reducible in 
polynomial time to the next one in the list, and GP is Turing reducible to WP 
in exponential time (see definitions in Section [4. 1[) . 

Among general facts on computational complexity of geodesies notice that 
if G has a polynomial growth, i.e., there is a polynomial p{n) such that for each 
n £ N cardinality of the ball Bn of radius n in the Cayley graph r(G, X) is at 
most p(n), then one can easily construct this ball Bn in polynomial time with an 
oracle for WP in G. If, in addition, such a group G has WP decidablc in polyno- 
mial time then all the problems above have polynomial time complexity with re- 
spect to any finite generating set of G (since the growth and WP stay polynomial 
for any finite set of generators). Now, by Gromov's theorem [25], groups of poly- 
nomial growth are virtually nilpotent, hence linear, so they have WP decidable 
in polynomial time. It follows that all Geodesic Problems are polynomial time 
decidable in groups of polynomial growth (finitely generated virtually nilpo- 
tent groups). On the other hand, there are many groups of exponential growth 
where GP is decidable in polynomial time, for example, hyperbolic groups [15] 
or metabelian Baumslag-Solitar group BS{l,n) = {a,t \ t~^at = a"), n > 2 
(see [13] and Section HT] for comments). 

In general, if WP in G is decidable in polynomial time then BGLP is in the 
class NP, i.e., it is decidablc in polynomial time by a non-deterministic Turing 
machine. It might happen though, that BGLP in a group G is as hard as any 
in the class NP, i.e., it is NP-complete. The simplest example of this type is 
due to Perry, who showed in [33] that BGLP is NP-complete in the metabelian 
group Z2wr(Z x Z) (the wreath product of Z2 and Z x Z). Correspondingly, the 
search problems GP and GLP are NP-hard, which means, precisely, that some 
(any) NP-complete problem is Turing reducible to them in polynomial time. 
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Our view-point on geodesies in free solvable groups is based on geometric 
ideas from the following two papers. In 1993 Droms, Lewin, and Servatius in- 
troduced a new geometric approach to study WP and GLP in groups of the 
type F/N' via paths in the Cayley graph of F/N [12]. In 2004 Vershik and Do- 
brynin studied algebraic structure of solvable groups, using homology of related 
Cayley graphs [SSj . This approach was outlined earlier in the papers [STJ [52] , 
where possible applications to random walks on metabelian groups have been 
discussed. In the papers jSTJ [52] (see also |53]) a new robust presentation of a 
free metabelian group Sr,2 was introduced as an extension of Z'' by the integer 
first homology group of the lattice Z'' (viewed as a one-dimensional complex) 
with a distinguished 2-cocycle. Similar presentations of other metabelian and 
solvable groups laid out foundations of a new approach to algorithmic problems 
in solvable groups. 

It seems these ideas are still underdeveloped in group-theoretic context, de- 
spite their obvious potential. Meanwhile, in semigroup theory similar geometric 
techniques have been widely used to deal with free objects in semidirect products 
of varieties. One can find an explicit exposition of these techniques in the pa- 
pers due to Almeida [1] and Almeida and Weil [2], while in |48l ?,[11[3] Auinger, 
Rhodes and Steinberg use similar machinery on a regular basis. Earlier, similar 
methods, though sometimes implicitly, were used in inverse semigroups theory, 
we refer here to papers [42 ] [37 ] [38 l [8] . 

In group theory most of the results in this area relied on various forms of the 
Magnus embedding and Fox derivatives. The role that the Magnus embeddings 
play in varieties of groups was clarified by Shmelkin [49]. In [39] Matthews 
proved that the conjugacy problem (CP) in free metabelian groups is decidable, 
and Kargapolov and Remeslennikov generalized this to free solvable groups Sr,d 
[52] . A few years later Remeslennikov and Sokolov described precisely the image 
of F/N' under the Magnus embedding and showed that CP is rcsidually finite in 
Sr,d [47]. Wc refer to a survey [46] on algorithmic problems in solvable groups. 

In Sections [2^ and [2?5l we study elements of groups of the type F/N' via flows 
on the Cayley graph F of F/N. It turns out the flow generated by a word w G F 
on the graph F directly corresponds to the Fox derivatives of w in the group 
ring ZF/N. This simple observation links together the techniques developed in 
group theory for the Magnus embeddings with the extensive geometric and the 
graph-theoretic machinery for flows. Indeed, the set of geometric circulations 
(flows where the Kirchhoff law holds for all vertices, including the source and 
the sink) form a group which is naturally isomorphic to the first homology group 
Hi (F, Z) of F. In this content the geometric circulations on F represent precisely 
the 1-cycles of F (viewed as 1-complex). The classical result in homology theory 
describes Hi{T, Z) as the abclianization of the fundamental group 7ri(F), which, 
in this case, is isomorphic to the free group N. Putting all these together one 
has another geometric proof of the Fox theorem, as well as the description of 
the kernel of the Magnus embedding. 

In Section [2771 we describe geodesies in groups F/N' as Euler tours in some 
finite subgraphs of F generated by the supports of the flows of the elements of 
F/N' on F. The description is geometric, explicit, and it gives a natural way 
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to compute the geodesic length of elements. In this part geometric ideas seem 
unavoidable. However, this simplicity becomes treacherous when one concerns 
the efficiency of computations. 

We prove that BGLP (relative to the standard basis) is NP-complete even 
in Sr,2- Consequently, the problems GP and GLP are NP-hard in Sr,2- To show 
this we construct a polynomial-time reduction of the Rectilinear Steiner Tree 
Problem (RSTP), which is NP-complete, to BGLP in Sr,2- The necessary infor- 
mation on RSTP is outlined in Section 14.31 and the proof of the main theorem 
is in Section [44l Notice, that in [12] GLP was claimed to be polynomial time 
dccidablc in arbitrary finitely generated free solvable groups, but the argument 
turned out to be fallacious. 

In the second half of the 20th century free solvable groups, as well as, solv- 
able wreath products of groups and finitely generated metabelian groups, were 
intensely studied, but mostly from the view-point of combinatorial groups the- 
ory. Now they stand at the heart of research in various areas of algebra. On the 
one hand, the rejuvenated interest to these groups stems from random walks 
on groups and, cohomology theory. For example, wreath products of abelian 
groups give exciting examples and counterexamples to several conjectures on the 
numerical characteristics of random walks. It seems, the main reasons that fa- 
cilitate research here come from some paradoxical properties of the groups itself: 
all these groups are amenable (as solvable group), but they have exponential 
growth and may have nontrivial Poisson boundary [31| . etc. These groups, con- 
trary to, say, free nilpotent groups, may have irreducible unitary representations 
with nontrivial cohomology. Some numerical characteristics of these groups are 
very intriguing, giving new exciting examples in the quantitative group theory. 
For example, metabelian "lamplighter" groups have intermediate growth of the 
drift, positive entropy, etc. These groups were intensively studied recently (see 
papers [31] [16] and the bibliography in the latter) . 

On the other hand, metabelian groups are currently at the focus of a very 
active research in geometric groups theory. In 1983 Gromov proposed a program 
for studying finitely generated groups as geometric objects |26j . One of the 
principal directions of this program is the classification of finitely generated 
groups up to quasi-isometry. It follows from Gromov's result on groups with 
polynomial growth [25] that a group quasi-isometric to a nilpotent group is 
virtually nilpotent. In the case of solvable groups the situation is much less 
known. Erschler shown in [13j that a group quasi-isometric to a solvable group 
may be not virtually solvable. Thus, the class of virtually solvable groups is not 
closed under quasi-isometry. On the other hand there are interesting classes of 
solvable non-polycyclic groups that are quasi-isometrically rigid, for example, 
solvable Baumslag-Solitar groups (Farb and Mosher [miH]). We refer to the 
papers [41] and [17] for some recent results in this area. 

It seems timely to try extend the results of this paper to the classes of 
solvable groups mentioned above. There are many interesting open questions 
concerning computational complexity of algorithmic problems in these classes 
of solvable groups, we discuss some of them in Section [5] 

All polynomial time algorithms presented in this work are implemented and 
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available at [TU] , 
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2 Preliminaries 

2.1 The Word Problem 

Let F = Ft — F{X) be a free group with a basis X = {xi, . . . , Xr}- A subset 
R C F defines a presentation P = {X \ R) oi a group G = F/N where N = 
ncl{R) is the normal closure of R in F. If R is finite (recursively enumerable) 
then the presentation is called finite (recursively enumerable). 

The Word Problem WP for P is termed decidable if the normal closure N is 
a decidable subset of F{X), i.e., there exists an algorithm A to decide whether a 
given word w G F{X) belongs to N or not. The time function : F{X) N 
of the algorithm A is defined as the number of steps required for A to halt 
on an input w E F(X). We say that the Word Problem for P is decidable in 
polynomial time if there exists a decision algorithm A, as above, and constants 
c,k&N such that 

TAiw) < c\w\'' 

for every w € F{X) (here \w\ is the length of the word w). In this case we say 
that the time complexity of WP for P is 0{n^). 



2.2 Free solvable groups and the Magnus embedding 

For a free group F = F{X) of rank r denote by i^^^^ = F' = [F, F] the derived 
subgroup of F, and by F^'^'^ ~ ^p{<i-i) ^ p^d-i)-^ _ ^j^g derived subgroup of 
F, d > 2. The quotient group Ar = Fr/F^ is a free abelian group of rank r, 
Mr = Fr/ Fr"^^ is & free metabelian group of rank and Sr^d = Fr/Fr'^'^ is a free 
solvable group of rank r and class d. In the sequel we usually identify the set X 
with its canonical images in A^jM^ and Sr^d- 

One of the most powerful approaches to study free solvable groups is via the, 
so-called, Magnus embedding. To explain we need to introduce some notation. 
Let G = F/N and ZG the group ring of G with integer coefficients. By /i : 
F —t G we denote the canonical factorization epimorphism, as well its linear 
extension to /i : ZF ZG. Let T be a free (left) ZG-module of rank r with a 
basis {ii, . . . , tr}- Then the set of matrices 



M{G)=\^ 1 j = U 1 j l.9€G,teT 

forms a group with respect to the matrix multiplication. It is easy to see that 
the group M(G) is a discrete wreath product M(G) = ArWrG of the free abelian 
group Ar and G. 
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In [33] Magnus showed that the homomorphism (p : F ^ M (G) defined by 




satisfies kercj) = N' . It foUows that (f> induces a mononiorphisni 

^ : F/N' ^ M{F/N), 

which is now called the Magnus embedding. 

The Magnus embedding allows one to solve WP in the group F/N' if WP 
in G = F/N is decidable. Indeed, given a word w G F{X) one can compute 
its image 4>{w) in M{G) (multiplying the images of the letters in w) and then, 
using a decision algorithm for WP in G, check if the resulting matrix 4){w) is 
the identity matrix or not. To estimate the complexity of such an algorithm 
notice first, that the coefficients from ZG that occur in the upper-right corner 
of the matrix have 0{\w\) summands. Secondly, to check whether or not 
an element ft. = miWi + . . . + mkVk G ZG, where nii ^ 'Z and Vi d G are given 
as words in the generators X from G, is trivial in ZG it requires about 0{P) 
comparisons of the type Vi = Vjl in G. This gives an estimate for the time 
function T of WP in F/N' via the time function T for WP in F/N: 

T'{n) = 0{rn^T{n)), 

where n = \w\. Since WP in Ar can be decided in linear time the estimate above 
shows that the complexity of WP in Mr is Oirn^). Moreover, induction on the 
solvability class d gives a polynomial estimate 0(r''-~'^n^''-~^) for WP in the 
free solvable group 5*,-,^. Thus, the Magnus embedding gives a straightforward 
polynomial time (in r and n) decision algorithm for WP in Sr,di but the degree 
of the polynomial grows with d. In particular, this algorithm is not polynomial 
as a uniform algorithm on the whole class of free solvable groups. 

2.3 Free Fox derivatives 

Let F = Fr{X) be a free group of rank r with a basis X = {xi, ...,Xr}- The 
trivial group homomorphism F ^ I extends to a ring homomorphism e : ZF 
Z. The kernel of e is called the fundamental ideal Ai? of ZF, it is a free (left) 
Zi^-module freely generated by elements xi — 1, . . . , — 1. 

In [^ni US HH [23] R. Fox introduced and gave a thorough account of the 
free differential calculus in the group ring ZF. Here we recall some notions and 
results referring to books [9] [5] [28] for details. 

A map D : ZF — > ZF is called a derivation if it satisfies the following 
conditions: 

(Dl) D{u + v) = D{u) + D{v); 

(D2) D{uv) = D{u)v^ + uD{v), where e is the ring homomorphism defined 
above. 
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For every Xi € X there is a unique derivation, the so-called, a free partial 
derivative d/dxi, such that dxj/dxi = Sij, where 6ij is the Kronecker's delta. 
It turned out that for every u S ZF 

u ~ u' — du/dxi{xi — 1) + ... + du/dxr{xr — 1). (1) 

Since Ai? is a free ZF-modulc the equality ([1]) gives another definition of the 
partial derivatives. 

Condition (D2) implies the following useful formulas, that allow one to com- 
pute easily partial derivatives of elements of "ZF 

dxj'^ /dxi^ -5ijX~^. (2) 

and, hence, for a word w = x^^ . . . x\'^ € F{X) one has 

dwidx, = E^ixJ^ . . . x^rX {dxl^ /dx,) = (3) 

— -^ii ■ ■■ -^ij^i -^11 ■ ■ ■ -^ij 

l<j<n, ij—i, Sj—l ^l^jl^^^ ^j— *j s:^— — 1 

The following result is one of the principle technical tools in this area, it 
follows easily from the Magnus embedding theorem, but in the current form it 
is due to Fox [20l EH [22l ES] . 

Theorem [Fox] Let N be a normal subgroup of F and n : F ^ F/N the 
canonical epimorphism. Then for every u € F the following equivalence holds: 

\/i {du/dx^Y = u G [N,N]. 



In particular, for A'' = F^"^^ the standard epimorphism fi : F ^ Sd = F/F^'^'^ 
gives rise to a ring homomorphism /i : ZF — > ZS'd such that 

F^'^+i) = {ueF\ {du/dxif = /or i = 1, . . . , r}. (4) 

Composition of d/dxi with ^ gives an mdwced partial derivative jdxi : TiF 
ZiSd, which we often denote again by d/dxi omitting fj, (when it is clear from 
the context). 

Partial derivatives d'^/dxi are useful when computing images under the Mag- 
nus embedding. Indeed, by induction on the length of w S F it is easy to show 
that the image of w, under the Magnus embedding (j> : F/N' AI{F/N), can 
be written as follows 

^ { ELi d^w/dxi ■ ti 
\0 1 

This shows that the faithfulness of the Magnus embedding is, in fact, equivalent 
to the Fox Theorem above. 



Word Problem in free solvable groups 



10 



2.4 Flows on F/N 

In this section we relate flow networks on the Caylcy graph of F/N to the 
elements of F/N'. 

Let X = {xi, . . . , Xr} be a finite alphabet. An X-labeled directed graph T (or 
X- digraph) is a pair of sets {V, E) where the set V is called the vertex set and 
the set E (IV xV X X is called the edge set. An element e = {vi,V2,x) G E 
designates an edge with the origin vi (also denoted by o(e)), the terminus V2 
(also denoted by t{e)), labeled by x. If for e £ £^ we have o(e) = t{e) then we 
say that e is a loop. The graph F can be finite or infinite. 

Example 2.1. The Caylcy graph r{G,X) of the group G = F/N is an X- 
digraph. 

Given an A"-digraph T, wc can make F into a directed graph labeled by 
the alphabet X^^ = X U X^^. Namely, for each edge e = {vi,V2,x) of F we 
introduce a formal inverse e^^ = {v2,vi,x~^). For the new edges we set 
(e~^)~^ = e. The new graph, endowed with this additional structure, is denoted 
by F. In fact in many instances wc abuse notation by disregarding the difference 
between F and F. 

Remark 2.2. If X is a generating set of G such that X n X~^ = then f (G, X) 
is the Caylcy graph F(G', A'^-'^) of G relative to the generating set X^^ . 

The edges of F inherited from F are called positively oriented or positive. The 
formal inverses of positive edges in F are called negatively oriented or negative. 
The edge set of f sphts in a disjoint union £'(f) = £'+(F) U E^ {T) of the sets 
of positive and negative edges. 

The use of F allows us to define the notion of a path in F. Namely, a path 
p in F is a sequence of edges p = ei, . . . , where each is an edge of F and 
the origin of each (for i > 1) is the terminus of e^-i. In this situation we say 
that the origin o{p) of p is o(ei) and the terminus t{p) is t{ek). The length \p\ 
of this path is set to be k. Also, such a path p has a naturally defined label 
v{p) = v{ei) . . . v{ek). Thus v{p) is a word in the alphabet Tj ~ X [JX~^ . Note 
that it is possible that v^p) contains subwords of the form aa^^ or a^^a for 
some a G X . If w is a vertex of F, we will consider the sequence p = v to be a 
path with o{p) = t{p) = v, \p\ = and v{p) = 1 (the empty word). 

In general, one can consider labels in an arbitrary inverse semigroups, the 
construction above applies to this case as well. In particular, we will consider 
directed graphs with labels in Z. We consider also digraphs with no labels at 
all (to unify terminology, wc view them sometimes, as labeled in the trivial 
semigroup {1}), the construction above still applies. 

Let F = (1/, E) be an X-digraph with two distinguished vertices s (called 
source) and t (called sink) from V . Recall that a flow (more precisely Z-flow) 
on F is a function f : E ^ 1 such that 

(F) for all I. e F - {s, t} the equality Y.oie)=. /(e) " Y.tie)=v /(e) = holds. 
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The number f*{v) = Y.o{e)=v fi^)-Y.t{e)=v fi^) called the net flow at v G V. 
The condition (F) is often referred to as the Kirchhojf law (see, for example, 
[3 HI]) or a conservation law [?]■ 

For the digraph F the definition above can be formulated in the following 
equivalent way, which is the standard one in flow networks: 

(Fl) /(e) = -/(e-i) for any e e E. 

(F2) Eo(e)=. /(e) - for all e 1/ - {s, t}. 

Here the net flow at v is equal to f*{v) — J2o{e)=v fi^)- 

Usually a flow network comes equipped with a capacity function c : E ^ N, 
in which case a flow / has to satisfy the capacity restriction 

(F3) /(e) < c(e) for aU e e E. 

In the sequel we do not make much use of the capacity function (it occurs in 
an obvious way), so in most cases we consider flows on graphs F satisfying the 
Kirchhoff law (F) (or, equivalently, on graphs F satisfying (Fl) and (F2)). 

A flow / is called a circulation if (F) holds for all vertices from V (including 
s and t). 

Example 2.3. Let F = T{G,X) be the Cayley graph of G = F/N relative 
to the generating set X. The constant function / : E{T) {1} defines a 
circulation on F, since for every vertex g 6 V{T) and every label x € X there is 
precisely one edge {gx~^,g) with label x incoming into g and precisely one edge 
(g,gx) with the label x leaving g. 

An important class of flows on F = F(G, X) comes from paths in F. A path 
p in F defines an integer-valued function TTp : E{T) — > Z, such that for an edge 
e 7r(e) is the algebraic sum (with respect to the orientation) of the number of 
times the path p traverses e, i.e., each traversal of e in the positive direction 
adds +1, and in the negative direction adds —1. It is obvious that ttp is a flow 
on F with the source o{p) (the initial vertex of p) and the sink t{p) (the terminal 
vertex of p). Notice that TTp satisfies also the following conditions: 

(F4) cither tt^ is a circulation (iff p is a closed path), or /*(s) ~ 1, f*{t) = —1. 

(F5) TTp has finite support, i.c, the set supp{tt) = {e ^ E \ 7r(e) ^ 0} is finite. 

We say that a flow tt on F is geometric if it satisfies conditions (F4) and (F5). 

It is easy to see that the set C(F) of all circulations on F forms an abelian 
group with respect to the operations (here /, 5 G C(F)): 

• (/ + .9)(e) = ,/(e)+.g(e). 



• i-f)ie) - -/(e). 
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Meanwhile, the set QC{T) of all geometric circulations is a subgroup of C(r). 
On the other hand, the sum f + g oi two geometric flows gives a geometric flow 
only if the sink of / is equal to the source of g, or either f oi g (or both) is a 
circulation. In fact, the set G^{r) of all geometric flows is a groupoid. 

Let n(r) be the fundamental groupoid of paths in F. Then the map a : 
n(r) G^(r) defined for p S n(r) by a{p) = tt^ is a morphism in the category 
of groupoids. i.e., the following holds (here p,q E n(r)): 

• = ■''"p + ^g: if PQ is defined in n(r). 

Now wc will show that every geometric fiow tt on F can be realized as a path 
flow TTp for a suitable path p. 

Lemma 2.4. Let -k be a geometric flow on T . Then there exists a path p in T 

such that TT = TTp. 

Proof. Let tt be a geometric flow on F with the source s and the sink t. Denote 
by F^ the subgraph of F generated by supp{7r) U {s, t}. Suppose Q is a subgraph 
of F such that A = FttUQ is a connected graph (every two vertices are connected 
by a path in F). Clearly, tt induces a flow on A. Now we construct another X- 
digraph A* by adding new edges to A in the following manner. For every edge 
e G E{A) with |7r(e)| > 1 we add extra |7r(e)| - 1 new edges e'^^\ . . . ,el''('=)l~^ 
from o(e) to t{e) if 7r(e) > and from t{e) to o(e), if 7r(e) < 0. We label 
the new edges by the same label if 7r(e) > 0, and by its inverse, otherwise. If 
7r(e) = then we add a new edge e~^ from t{e) to o(e) with the inverse label. 
In the case |7r(e)| = 1 we do not add any new edges. Notice that every vertex in 
V{A*) — {s, t} has even directed degree (the number of incoming edges is equal 
to the number of outgoing edges). There are two cases to consider. 

Case 1). Suppose tt is a circulation. Then every vertex in A* has even 
directed degree. Therefore, the digraph A* has an Euler tour p* , i.e., a closed 
path at s that traverses every edge in A* precisely once. Let : A* — > A be the 
morphism of X-digraphs that maps all the new edges e'^\ . . . , gl'^^'^'l^^ to their 
original edge e. Clearly, the image p = 4>{j)*) is a path in A such that Hp = tt. 

Case 2). Suppose tt is not a circulation. Let g be a path in A* from s to t. 
Then tt' = tt — tt^ is a circulation. Hence by Case 1) there exists a path p in F 
such that tt' — Tip. Therefore, Ttpq = -Kp + Hq — tt' + -Hq — tt, as required. 

□ 

2.5 Geometric interpretation of Fox derivatives 

In this section we give a geometric interpretation of Fox derivatives. 

Let G = F/N, p : F F/N the canonical epimorphism, and F = F(G', X) 
the Cayley graph of G with respect to the generating set X . A word w S F{X) 
determines a unique path p^ in F labeled by w which starts at 1 (the vertex 
corresponding to the identity of G) . As we mentioned in Section 12.41 the path 
Pw defines a geometric flow 7rp^„ on F which wc denote by tt^,. 
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Lemma 2.5. Let w € F = F{X). Then for any g € F/N and x £ X the value 
o/tTu, on the edge e ~ {g,gx) is equal to the coefficient in front of g in the Fox 
derivative (dw/dx)'^ € IG, i.e., 

(dwjdxY ^ ^ TT.u,{g,gx)g 
geG.xex 

Proof. Follows by induction on the length of w from formulas jS]). □ 

Figure [T] is an example for F — F{{xi, X2}) and N ~ F' . Non-zero values of 
TTw are shown as weights on the edges (zero weights are omitted). 



-10 12 3 



Figure 1: The values of Tr^, for w = a;2a;ia;2a;ia;2a^i ^x^ ^ on (xi, a;2)-grid. In 
this case dw/dxi ~ —\ + X2 — x\x\^x\x\ and dwjdxi = 1 — xi + x\x\ — x\x\. 

The following theorem has been proven in [121 ES] using a homological ar- 
gument similar to the one in Proposition l2.8l Here we give a short independent 
proof based on the Fox theorem. 

Theorem 2.6. \VA\5^ Let N he a normal subgroup of F and an equivalence 
relation on F defined by 

Then F/N' ^F/ r^N- 

Proof. Let u,v £ F. Suppose m = t; in F/N'. Then uv^^ G N' hence by 
Fox theorem d^{uv~^)/dx = for every x & X (here by /dx we denote the 
canonical image of d/dx in the group ring 7L(F /Nf). Observe that 

) = ^. (5) 



Hence in Z{F/N) 
so, by Lemma 12.51 7r„ = tt^,, as claimed 



dx dx dx ' 



Word Problem in free solvable groups 



14 



To show the converse, notice first that 7r„ = tTu implies that u = v in F/N. 
Indeed, it can be seen from the definition of tt but also follows from ([T]) and 
Lemma 12.51 since in this case 

K - f'") -{u- vY = V —(u -v)-(x-l) = 



x<^X 



which implies = . Now by ([5]) 



dx dx dx ' 

and, hence, by Fox theorem uv^^ G N' . 

□ 



Remark 2.7. Theorem 12.61 relates the algebraic and geometric points of view 
on derivatives. One can prove this theorem (see Section [2. 6p using a pure topo- 
logical argument, then the Fox theorem, as well as, the Magnus embedding, 
come along as easy corollaries. 



2.6 Geometric circulations and the first homology group 
of r 

We describe here geometric circulations on F = T{G,X) in a pure topological 
manner. For all required notions and results on homology of simplicial com- 
plexes we refer to [30] or |50j . 

In this section we view F as an infinite 1-complcx. 

Proposition 2.8. Let G = F/N, F = F(G,X), and a : 7ri(F) QC{T) be a 
map defined by (t{p) = TTp for p G 7ri(F) . Then: 

• cr is an epimorphism of groups; 

• every geometric circulation on F defines a 1-cycle on F; 

• Hi(F,Z) ~gC(F); 

• 7ri(F) ~ N and kercr = N' . 

Proof. It was mentioned already in Section 12.41 that the map p — > TTp is a mor- 
phism from the fundamental groupoid 11 (F) of paths in F into the groupoid of 
geometric fiows QT(T). Hence, the restriction of this map onto the fundamental 
group 7ri(F) of F gives a homomorphism of groups. We have seen in Lemma 1 2. 41 
that a is onto. This proves the first statement. 

To see 2) observe first that a geometric circulation / : E{r) Z, viewed 
as a formal sum X]ee-E(r) /(^)^' gives precisely a 1-chain in F (see, for example, 
[30]). Moreover, by definition, the net fiow /*(«) at the vertex v G V(r) is the 
coefficient in front of v in the boundary df of /. Therefore, df = 0, so / is a 
1-cycle. 
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3) follows easily from the 2). Indeed, T is an 1-complex, so there are no 
non-trivial 1-boundaries in T. In this event i?i(r, Z) is isomorphic to the group 
QC{T) of 1-cycles, as claimed. 

4) It is a classical result that the kernel of a : 7ri(r) — > Hi{T,'Z) is equal to 
the derived subgroup of 7ri(r) (see [30]). To prove 4) it suffices to notice that 
7ri(r) ~ N, which is easy. □ 

Remark 2.9. Proposition 12.81 gives a simple geometric proof of Theorem 12.61 
that relates the algebraic and geometric points of view on derivatives. Now one 
can derive the Fox theorem from the geometric argument above and then obtain 
the description of the kernel of the Magnus embedding, as an easy corollary. 

2.7 Geodesies in F/N' 

Let G = F/N and ^ : F{X) G the canonical epimorphism. In this section 
we describe geodesies of elements of the group H = F/N' relative to the set of 
generators X^^. 

It is convenient to view the free group F = F{X) as the set of all freely 
reduced words in the alphabet X^^ = XUX^^ with the obvious multiplication. 

To describe geodesies in H (relative to X) of a given word w G F we need a 
construction from Lemma 12.41 Recall that is the path in the Cayley graph 
r = r(G, X) from 1 to w'^ with the label w and tt^ the induced geometric flow 
on r with the source 1 and the sink w'^. By Fy, we denote the subgraph of 
F generated by supp{tt^) U {1, w'^}. Suppose Q is a finite subgraph of F such 
that A = Fu, U Q is a connected graph and Q has the least number of edges 
among all such subgraphs. It follows from minimality of Q that every connected 
component of Q is a tree. Moreover, if in the graph A = F^j U Q one collapses 
every connected component F^, to a point, then the resulting graph is a tree. 
We refer to Q as a minimal forest for w. In general, there could be several 
minimal forests for w. 

Similarly as in the proof of Lemma 12.41 we construct a finite X-digraph A* 
by rephcating every edge e G E{A) with ||7rt„(e)| — 1| new edges in such a way 
that every vertex in V{A*) — {1, w'^} has even directed degree and the map, that 
sends every replica of an edge e back to e (or e~^ depending on the orientation) 
is a morphism of X-labeled digraphs : A* A. There are two cases. 

• Case I. Suppose that p^ is a closed path in F, i.e., w G N. In this case 
every vertex in A* has even degree, so A* has a Euler tour, say Pq at 
1. Denote by pq the image 4>{p*q) of p*Q under (j). It follows from the 
construction (see Lemma [2^ that pq is a closed path at 1 in F such that 
T^w = T^pQ- Therefore, if Wp^ G F is the label of pq then by Theorem 12.61 
w = in H . Moreover, since pq is a Euler tour in A* its length, hence 
the length of Wp^ is equal to 

\pq\= E Me)\+2\E{Q)\- (6) 
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• Case II. Suppose that is a not a closed path in F, i.e., w ^ N. By 
induction on \w\ it is easy to show that the vertices 1 and w'^ belong 
to the same connected component of Fj^. Again, there exists an Euler's 
tour Pq in the graph A* which starts at the source and ends at the sink. 
Clearly, TTpg satisfies the equality If w is the label of the path tt^^ 
then 7r„ = tt^^ ~ 1:^ and w is a geodesic word for w. 

Now, with the construction in place, we are ready to characterize geodesies 
in H of elements from N . 

Theorem 2.10. Let H = F/N' and w £ F. Then the following holds: 

• if Q is a minimal forest for w then Wp^ is a geodesic for w and 

lx{w) = 

• every geodesic word for w is equal (in F) to a word Wp^ for a suitable 
minimal forest Q and an Euler path p*Q . 

Proof. Let u G F be a geodesic word for w^^ in H. Observe, that A = supp{'!Tu)U 
Pu is a connected subgraph of F and 

\Pu\> ^ 7r„(e) + 2|£'(A - supp(7r„))|. 

Now, by Theorem 12.61 the equality u = w in H implies tt^ = t^w In particular, 
suppiiTu) = supp{'Kyj). Hence 

\Pu\> T^u,{e) + 'i\E{A~ supp{^.^))\. (7) 

e^supp{p^, ) 

Since u is geodesic for if the number |A — supp{T:w)\ is minimal possible, so 
Q = A — supp{T:yj) is a minimal forest for w. In fact, the equation [7] shows that 
the converse is also true. This proves the theorem. □ 

The discussion above shows that CP is easy in F/N' provided one can solve 
the following problem efficiently. 

Minimal Forest Problem (MFP) Given a finite set of connected finite sub- 
graphs Fi, . . . , Fs in F find a finite subgraph Q of F such that Fi U . . . U F^ U Q 
is connected and Q has a minimal possible number of edges. 

Proposition 2.11. CP in F/N' (relative to X) is linear time reducible to MFP 
for T{F/N,X). 

Proof. Follows from the discussion above. Indeed, given a word w £ F one can in 
linear time compute the flow tt^, and find the connected components Fi, . . . , F^ 
of supp{-K^) U {o(p,„),t(p^)} in F = V{F/N,X). Then, solving MFP for these 
components, one gets the subgraph Q which makes the graph Fi U . . . U F^ U Q 
connected. Obviously, it takes linear time to find a Euler path in the graph A, 
hence to find a geodesic for w. □ 
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3 The Word Problem in free solvable groups 

In this section we present fast algorithms to compute Fox derivatives of elements 
of a free group F in the group ring liSr.d of a free solvable group Sr.d- As an im- 
mediate application, we obtain a decision algorithm for WP in a free metabelian 
group Mr with time complexity 0(rn log2 n) and a decision algorithm for WP in 
Sr,d,d > 3, with time complexity 0{rdn^). These are significant improvements 
in comparison with the known decision algorithms discussed in Introduction 
and Section 12.21 As another application we get a fast algorithm to compute 
images of elements from Sr.d under the Magnus embedding, which opens up 
an opportunity to use efficiently the classical techniques developed for wreath 
products. 

3.1 The Word Problem in free metabelian groups 

In this section we compute Fox derivatives of elements of F in the group ring 
Z(F/F'). Then we apply this to WP in free metabelian groups. 

Let X = {xi, . . . ,a;,.}, F = Fr = F{X), M = M- = F/F'^^\ A = Ar = 
F/F', and fi : F ^ A the canonical epimorphism. All Fox derivatives in this 
section are computed in the ring ZA. 

hetweF. Then 



One can encode all the derivatives d^w/dxi in one mapping Mw : A x 
{1, . . . , r} — i- Z, where Mw{a, i) = ma.i- Let 



and Sw the restriction of A/„, onto supp{Mm) . To compute Fox derivatives of w 
we construct a sequence of finite maps 5*0 = 0, 5*1, . . . , S'„ = S"™, as we read w. 
On each step k we either extend the domain Dom(Sk) of Sk or change the value 
of Sk on some element from Dom{Sk)- To do this we need a data structure 
which allows one to do the following operations efhciently: 

• for a given {a,i) determine if {a,i) G Dom(Sk) or not; 

• add (a, i) to Dom{Sk) if (a, i) ^ Dom(Sk) and define Sk{a, i) ~ q for some 



• change the value of Sk on (a, i) if {a,i) G Dom{Sk)- 

Every element a € A can be written uniquely in the form a = xj^^"'' . . . xi''^°'\ 
where (5,; (a) G Z, so one may use the tuple of coordinates (5(a) = 
{6i{a), . . . , 6r{a)) to represents a. It follows from the formula ([3]) for par- 
tial derivatives that for every (a,z) G Dom{Sw) the components Sj{a) of (5(a) 
satisfy the inequality |(5j(a)| < as well as the values of S^ and, hence, 




supp{Mw) = {(a,i) I Mru{a,i) ^ 0}, 



g G Z; 
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\Sw{a,i)\ < \w\. Therefore, it takes [logjd^l + 1)] bits to encode one coordi- 
nate Sj{a) (one extra bit to encode the sign + or rflogjdw] + 1)] bits to 
encode 6{a), and at most r[log2(|w| + 1)] + [log2 r] bits to encode (a,i). We 
denote the binary word encoding {a,i) by {a,i)*. 

Thus every function Sk can be uniquely represented by a directed {0, 1}- 
labeled binary tree with the root e and with leaves labeled by integers such 
that (a, i) G Dom{Sk) if and only if there exists a path in Tk from the root 
e to a leaf labeled by the code of (a,i) and such that the leaf of this path is 
labeled precisely by the integer Sk{a,i). Notice, that the height of Tk is equal 
to r[log2(|w| + 1)] + [logjr]. Such a tree is visualized schematically in Figure 

m 




r (log2 \w\+l) +log2r 

Figure 2: The tree Tk representing the function Sk- 

Remark 3.1. It is clear that one can perform every operation mentioned above 
on this data structure in at most r[log2(|w| + 1)] + [log2 r] elementary steps. 

Wc use this data structure to design the following algorithm for computing 

Algorithm 3.2. (Computing Fox derivatives in F/F') 

Input, r G N and w = x^^ . . . x^^^ e F{X), where ij S {1, . . . , r} and Ej = ±1. 

Output. Sw 

Computations. 

A. Set = and 6{a) = (0, . . . , 0) £ Z\ 

B. For j = 1, . . . ,n do: 

(1) if = 1 then 

* check if there is a path (from the root to a leaf) labeled by 
((5(a), ij)* in S; 

* if such a path does not exist in S then create it, add it to S, and 
put 1 as the corresponding value at the new leaf; 

* if such a path exists in S then add 1 to the value at its leaf. 
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(2) if Ej = -1 then 

* add Ej to the Zjth coordinate of <5(a); 

* check if there is a path (from the root to a leaf) labeled by 
((5(a), ij)* in S; 

* if such a path docs not exist in S then create it, add it to S, and 
put —1 at the corresponding leaf; 

* if such a path exists in 5* then subtract 1 from the value at its 
leaf. 

C. Output Sn- 

Theorem 3.3. Given r E N and w £ F Algorithm \3.S\ computes all partial 
derivatives of w in (the mapping S^) in time 0{r\w\ log2 \w\). 

Proof. Using the formula ^ for partial derivatives it is easy to check that 
given w £ F Algorithm 13. 2[ indeed, computes the mapping Sw. To verify the 
complexity estimates observe, first that there are iterations in the part [B.] of 
Algorithm l3.2l Each such iteration requires 0(r [log2(|iu| + 1)] ) elementary steps 
(see Remark 13. ip . so altogether one has 0{r\w\log2\w\) as the time complexity 
estimate for Algorithm 13.21 as claimed. 

□ 

Algorithm 3.4. (Word Problem in free metabelian groups) 

Input, r eN and w e F. 

Output. True if w represents the identity in Mr and False otherwise. 
Computations. 

• Apply Algorithm 13.21 to compute Sw 

• Check, looking at the values assigned to leaves of S^, if all Fox derivatives 
dw/dxi, i = 1, . . . , r are equal to or not. 

• If all the derivatives are then output True. If there is a non-zero deriva- 
tive then output False. 

Theorem 3.5. Algorithm \3.4\ solves the Word problem in a free metabelian 
group Mr in time 0{r\w \ log2 1^1). 

Proof. Follows from Theorem 13.31 and the Fox Theorem (see Section 12. 3p . □ 

3.2 The Word problem in free solvable groups 

In this section we present an algorithm to compute all Fox derivatives of a word 
w € F in the group ring ZS'^.d-i, d > 2 in time 0{rd\w\^). This gives a decision 
algorithm for WP in 5*^,^ within time complexity 0{rd\n\^). 

Let X = {xi,...,x'r}, F = Fr = F{X), S = = F/F^'^\ d > 3, and 
fi : F Sr,d-i the canonical epimorphism. All Fox derivatives in this section 
are computed in the ring ZSr.d-i- 
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In Section 13.11 we used a unique representation of elements a G Ar hy their 
coordinate vectors 6{a) to compute Fox derivatives in nearly linear time. Now 
we do not have such normal forms of elements of Sr,d, d > 2, so our computations 
are slightly different, however, the general strategy is quite similar. To speed 
up computations wc use some data structures based on efficient partitioning 
techniques. 

In general, let G be a group generated by X and D a finite subset of F{X). 
A G-partition of D is a partition of D into a union of disjoint non-empty subsets 
Di such that for any u,v G D, u ~ v in G if and only if they belong to some 
subset Di. Clearly, the G-partition of D is unique. Observe that if a group 
-ff is a quotient of G then the G-partition of D is the same or finer than the 
if -partition of D. 

If the set D is ordered, say D — {wq, . . . , w„} then the G-partition of D 
can be represented by a function P : {0, . . . , 71} — > {0, . . . , n} where P{j) = i if 
and only if Wi = Wj in G and i is the smallest with such property. Given the 
G-partition P of D one can arrange his data in such a way that it takes linear 
time (in the size of j) to compute P{j). In particular, given i,j one can check 
in linear time if Wi ~ Wj in G. Also, for a given word w Cz D one can find in 
linear time an index i such that w — Wi. These are the main two subroutines 
concerning partitions of D. 

Let w = . . . x^^ , where ij G {1, . . . , r} and ej ~ ±1. Put 

D.^ = {e, xll, xllxZ, xll...xl:]^F{Xr). (8) 

We order the set as follows iuq = e, . . . , = w. Now, to check whether 
or not the derivative dw/dxi is trivial in "LSr^d-i one has to determine which 
pairs {wi,Wj) of elements from represent the same element in Sr.d-i a-nd 
then cancel out Wi with Wj in dw/dxi if they have the opposite signs. 

The goal of Algorithm 13. 101 below is to compute S'^^d-partition for D.^. This 
is performed in a sequence of iterations. The algorithm starts out by computing 
the A^-partition of Z?^,. On the second iteration the algorithm computes the 
M^-partition of On the third step it computes the 5'r,3-partition of D^]- It 
continues this way until it computes the ^'r-^d-partition of D^]. 

To explain how the algorithm works assume that the S'r,d-i-partition of 
is given by the partition function Pd^i described above. Notice, that the 5*1.,^- 
partition of D^^ is the same or finer than the S'^.d-i-partition of D^j, since Sr.d-i 
is a quotient of Sr,d- This shows that to construct S'r,(i-partition Pd of one 
has only to compare elements from which are equal in Sr,d-i- Suppose 
that Ws,Wt € -Dto, s < t and Ws = wt in 5*^,^-1. To check if Wg ~ wt in Sr,d 
we compare all their Fox derivatives in Z5r.d-i, so for every k = 1, . . . ,r we 
compute the following differences: 

dws/dxk - dwt/dxk = 

/ J <■! '■J — 1 / -/ tl 

l<j<s, ij—k, l^i^-^i ij—k. ej—~l 
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l<j<t, ij—k, £j— 1 1^^^^: ij—k, ej— — l 



s+l<j<t, ij—k, ej—1 '^ + 1^:?^^! ij—k. Sj — — ! 



E E 

.s + l<j<i, ij—k. Sj — l s+l<j<i, ij—k, ej—~l 

Clearly, given w and s, as above one can compute the formal expression ([9]) in 
time 0(|i(;|). To check if viewed as an element in ZSr,d-i, is equal to 
it suffices to represent it in the standard group ring form X^geS d i (where 
m G Z) and verify if all coefficients in this representation are zeros. Now we 
describe a particular procedure, termed Collecting Similar Terms Algorithm, 
which gives the standard group ring form for Q. Given ^ one can compute 
in time 0(|w|) the following sum 

E ^pu-i) + E ^po")' (10) 

s+l<j/<f, ij—k, Cj—l s-|-l<j/<#, ij—k, Cj— — ! 

Observe now, that two summands Wp and Wq in pop are equal in Sr.d-i if a-nd 
only if p = q. It is easy to see that it takes time 0(|w|) to collect similar terms 
in (|10p . i.e., to compute the coefficients in the standard group ring presentation 
of (Uni). 

It follows that dws/dxk = dwt/dxk in ZS'r,^-! if and only if all the co- 
efficients in the standard group ring form of pop are equal to 0. The argu- 
ment above shows that one can check whether or not dwg/dxk = dwt/dxk 
in liSr.d-i in time 0(|w|). Since we need to compare all partial derivatives 
d/dxk, fc = 1, . . . , r, of the elements Wg and wt it takes altogether 0(r|i(;|) time 
to verify if Ws = Wt in Sr.d-i- 

The routine above allows one to construct effectively the S'^.d-partition Pd 
of when given the iS'r,d_i-partition Pd-i- A more formal description of the 
algorithm is given below. 

Algorithm 3.6. (Computing the S'^.d-partition of 

Input. Positive integers r > 2, > 2, and a word w £ F{X). 
Output. The S'r,d-partition function Pd for the set Dyj- 

Initialization. Compute the set Dyj and form the initial (trivial) F-partition 

Po of Du,, so Po{i) = i for every z e {0, . . . , n}. 

Computations. 

A. Compute A-partition Pi of 

B. For c = 2, . . . , d do: 

1) For each Q < s < t < n such that Wg ~ Wt in Sr,c-i check whether 
or not Ws = Wt in Sr,c- 

2) Form the S'^.c-partition Pd of 
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C. Output Pd. 

Lemma 3.7. Given integers r,d > 2 and w £ F Algorithm Iff. 61 computes the 
Sr,d-partition (the function Pd) of D^t, in time 0{dr\w\^). 

Proof. Algorithm 13.61 makes precisely d iterations c = 1, . . . , d by consequently 
computing the 5r,c-partitions of Du,. After the S^.c-i-partition of D^j is com- 
puted, the algorithm computes the S'rx-partition of by comparing elements 
Ws, Wt G in Sr,c- It requires at most + l)/2 such checks, and, as was 

explained above, every such check can be done in time 0(r|w|). Altogether one 
needs 0(r|wp) steps to construct the function Pc on the iteration c. Since the 
algorithm makes altogether d iterations it takes it 0{dr\w\^) time to produce 
Pd, as claimed. 

□ 

Now we are in a position to show two applications of Algorithm 13.61 The 
first one concerns with computing Fox derivatives in 'LSr^d- 

Algorithm 3.8. (Computing Fox derivatives) 

Input. Positive integers r > 2, d > 2, a word w G F{X), and a number 

fce{i,...,r}. 

Output. The standard group ring presentation of the Fox derivative dw/dxk 

in ZSr,d- 

Computations. 

A. Compute, using formals the Fox derivative dw/dxk in ZF. 

B. Compute, using Algorithm 13.61 the S'^.^-partition of Dy^. 

C. Compute, using the Collecting Similar Terms Algorithm, the standard 
group ring form of dw/dxk in 'ZSr,d 

D. Output dw/dxk computed in [C.]. 

Lemma 3.9. Given integers r,d > 2, a word w £ F, and a number k g 
{1, . . . , r} Algorithm \3.8\ computes the standard group ring presentation of the 
Fox derivative dw/dxk in "LSr^d in time 0(dr\w\'^). 

Proof. Follows from Lemma 13.71 □ 

The second application of Algorithm 13.61 is to WP in Sr,d- 
Algorithm 3.10. (WP in Sr^) 

Input. Positive integers r > 2, d > 2, and a word w £ F{X). 
Output. True if ui = 1 in Sr,d and False otherwise. 
Computations. 

A. Compute the set Duj. 

B. Using Algorithm 13.61 compute 5'r,d-partition Pd of D^. 
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C. If Pdi^) = Pdin), i.e., 1 = ui in Sr4, then output True. Otherwise output 
False. 

Theorem 3.11. Algorithm \3.1(A solves the Word Problem in a free solvable 
group Sr,d 'in time 0{dr\w\^). 

Proof. Follows immediately from Lemma 13.71 □ 

4 Geodesies in free metabelian groups 

In this section we discuss the computational hardness of different variations of 
Geodesic problems and prove the main result about NP-completeness of BGLP 
in free metabelian groups. 

4.1 Algorithmic problems with geodesies in groups 

Let G be a group with a finite set of generators X = {xi, . . . , x^} and \x : 
F{X) — > G the canonical epimorphism. In this section we view the free group 
F{X) as the set of all freely reduced words in the alphabet X^-^ = X U X~^ 
with the obvious multiplication. 

For a word w in the alphabet X^-^ by \w\ we denote the length of w. The 
geodesic length of an element g € G relative to X, denoted by lx{g), is the 
length of a shortest word w G F{X) representing g, i.e., 

lx{g) = min{\w\ I w e F{X), = g}. 

To simplify notation we write, sometimes, Ixi'w) instead of Ixiw'^). A word 
w G F{X) is called geodesic in G relative to X, if \w\ = Ixiw). 

We are interested here in the following algorithmic search problem in a given 
group G described as above. 

The Geodesic Problem (GP): 

Given a word w e FiX) find a geodesic (in G) word w G F{X) such that 
w^^ = uj'^. 

One can consider the following variation of GP. 

The Geodesic Length Problem (GLP): 

Given a word w G F{X) find lx{w). 

Though GLP seems easier than GP (since a solution to GP gives, in linear 
time, a solution to GLP), in practice, to solve GP one usually solves GP first, 
and only then computes the geodesic length. 

As customary in complexity theory one can modify the search problem GLP 
to get the corresponding bounded decision problem: 

The Bounded Geodesic Length Problem (BGLP): 
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Let G be a group with a finite generating set X. Given a word w S F{X) and 
a natural number k determine if lx{w) < k. 

It is instructive to compare the algorithmic "hardness" of the problems above 
and the Word Problem (WP). Clearly, if one of them is decidable then all of 
them are decidable. To see the difference we need to recall a few definitions. 
Let A and B be algorithmic problems with inputs sets Ia and Ib- Then A is 
termed Turing reducible to B in polynomial time, if there exists an algorithm 
A with an oracle for B (which can be viewed as a "subroutine" of A that for 
a given input e G in one step returns the answer for B on e) that solves A 
in polynomial time. Similarly, one can define Turing reducibility in exponential 
time. In these cases we write A :<T.p B or, correspondingly, A diT,exp B. 

Again, it is not hard to see WP diT,p BGLP ^t,p GLP ^t,p GP. More- 
over, since (by brute force algorithm) GP d:T,exp WP it follows that all these 
problems are Turing reducible to each other in exponential time. Moreover, if 
G has polynomial growth, i.e., there is a polynomial p{n) such that for each 
nN cardinality of the ball i?„ of radius n in the Cayley graph r(G, X) is at 
most p(n), then one can easily construct this ball Bn in polynomial time with 
an oracle for the WP in G (see, for example, [12] for details). It follows that 
if a group with polynomial growth has WP decidable in polynomial time then 
all the problems above have polynomial time complexity with respect to any 
finite generating set (since the growth and WP stay polynomial for any finite 
set of generators). Observe now, that by Gromov's theorem [25| finitely gener- 
ated groups of polynomial growth are virtually nilpotent. It is also known that 
the latter have WP decidable in polynomial time (nilpotent finitely generated 
groups are linear). These two facts together imply that the Geodesic Problem 
is polynomial time decidable in finitely generated virtually nilpotent groups. 

On the other hand, there are many groups of exponential growth where GP 
is decidable in polynomial time, for example, hyperbolic groups [,15j . Among 
metabelian groups the Baumslag-Solitar group i?S'(l,2) = {a,t \ t^^at = a?) 
has exponential growth (it is solvable but not polycyclic - and the claim follows 
from the Milnor theorem [3D]) and GP in _B5'(1,2) is decidable in polynomial 
time (see [T3]). 

In general, if WP in G is polynomially decidable then BGLP is in the 
class NP, i.e., it is decidable in polynomial time by a non-deterministic Turing 
machine. Indeed, if Ixiw) < k then there is a word u e F{X) of length at most 
k which is equal to w in G - this m is a "witness" of polynomial size which allows 
one to verify in polynomial time that lx{w) £ k (just checking that u = w vnG). 
In this case GLP is Turing reducible in polynomial time to an NP problem, but 
we cannot claim the same for GP. Observe, that BGLP is in NP for any finitely 
generated metabelian group, since they have WP decidable in polynomial time 
(see Introduction). 

It might happen though, that WP in a group G is polynomial time decidable, 
but BGLP in G is as hard as any problem in the class NP, i.e., it is NP-complete. 
Recall (in the notation above) that a decision problem B is NP-complete if it is in 
NP and for any decision problem A from NP there is a computable in polynomial 
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time function f : Ia ^ Ib (Karp reduction, or a polynomial reduction), such 
that A is true on x € Ia if and only if B is true on f{x). The simplest example 
of this type is due to Perry, who showed in [32] that BGLP is NP-complete in 
the metabelian group l^wril x Z) (the wreath product of Z2 and Z x Z). In 
this event, the search problems GP and GLP are called NP-hard, this means 
precisely that some (any) NP-complete problem is a Turing reducible to them 
in polynomial time. 

It would be very interesting to classify finitely generated metabelian groups 
with respect to computational hardness of their GP or GLP problems. In the 
next section we clarify the situation with free metabelian groups. Some remain- 
ing open problems are discussed in Section [S) 

It was claimed in [12] that in free solvable groups of finite rank GLP is dccid- 
able in polynomial time. Unfortunately, in this particular case their argument 
is fallacious. Our main result of this section is the following theorem. 

Theorem 4.1. (Main Theorem) Let Mr be a free metabelian group Mr of 
finite rank r > 2. Then BGLP in Mr (relative to the standard basis) is NP- 
complete. 

Proof. The proof of this result consists of two parts. Firstly, in Section 14.21 
(Corollary [53]) wc show that it suffices to prove that BGLP is NP-complete in 
M2. Secondly, in Section (Theorem I4.11|) wc give a proof that BGLP is, 
indeed, NP-complete in M2. □ 

This immediately implies the following results. 

Corollary 4.2. The search problems GP and GLP are NP-hard in non-abelian 
Mr (relative to the standard basis). 

To prove the Main Theorem wc reduce the problem to the case r — 2 and 
then show that BGLP in M2 is NP-complete. To sec the latter wc construct a 
polynomial reduction of the Rectilinear Steiner Tree Problem to BGLP in Mr. 

4.2 Reduction to M2 

Let V be a variety of groups. For groups A, B E V we denote hy A *\; B the 
free product of A and B relative to V. In particular, if A = {X \ R), and 
B = {Y \ S) are presentations of A and B in V then A*v B ^ {xUY \ RU S) is 
a presentation oi A*v B in V. As usual, A*\> B satisfies the canonical universal 
property: any two homomorphism A C, B —^ C into a group C £ V extends 
to a unique homomorphism A*y B ^ C (we refer to [33 for details). It follows 
that if Fv{X U F) is a free group in V with basis X UY then F^iX UY) = 
FviX)^vFv{Y). 

The following lemma claims that free V-factors of a group G are isometrically 
embedded into G. 

Lemma 4.3. Let A, B £ V with finite generating sets X and Y. Then in the 
group A*\;B no geodesic word (relative to XUY) for an element from A contains 
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a letter from Y . In particular, for any word w € F{X) its geodesic length in A 
(relative to X) is equal to the geodesic length in A*xi B (relative to X [JY ). 

Proof. Let w £ F{X) be a geodesic word in A relative to X. Suppose that 
u G F{X U F) is a geodesic word vci G = A B (relative to X \J Y) that 
defines the same element as w. The identical map A A and the trivial map 
B — > 1 give rise to a homomorphism (p : G ^ A. This cf), when applied to u, 
just "erases" all letters from Y . It follows that if u contains a letter from Y 
then \u'^\ < \u\ < \w\ - contradiction with the assumption that w is geodesic in 
A relative to X (since w — u'^ in A). □ 

Corollary 4.4. In the notation above, each of the problems GP, GLP, BGLP 
in A (relative to X) is polynomial time reducible to the problem of the same 
type in A*\> B (relative to X UY). 

Notice, that Mn+m = Mn *A2 M„i, where A2 is the variety of all metabelian 
groups. Since WP is in P for groups from A4 Corollary implies the following 
result. 

Corollary 4.5. If BGLP is NP-complete in M2 then it is NP-complete in 
Mr, r > 2, relative to the standard bases. 

Remark 4.6. Corollary 14.51 easily generalizes to free groups in an arbitrary 
variety, provided they have WP decidablc in polynomial time. 

4.3 Rectilinear Steiner Tree Problem 

The Steiner Tree Problem (STP), which was originally introduced by Gauss, is 
one of the initial twenty one NP-complete problems that appeared in the Karp's 
list [33] . We need the following rectilinear variation of STP. 

Let be the Euclidean plane and F the integer grid canonically embedded 
into ffi^ (all vertices from together with all the horizontal and vertical lines 
connecting them). If A is a finite subset of then a rectilinear Steiner tree 
(RST) for A is a subgraph T of F such that AUT is connected in F. By s{T) 
(size of T) we denote the number of edges in T. An RST for A is optimal if it 
has a smallest possible size among all RST for A, we denote such RST by Ta- 
Observe, that a given A may have several different optimal RST, but their size 
is the same, we denote it by s{A). 

Notice that, in general, Ta for A is not a spanning tree for A in T, it may 
contain some vertices from which are not in A (so-called, Steiner points). 
The Rectilinear Steiner Tree Problem (RSTP) asks for a given finite A C 
and fe G N decide if there exists some Ta for A with s{Ta) < k. It is known 
that RSTP is NP-complete [M]. 

4.4 NP-completeness of BGLP in M2 

Now we construct a polynomial reduction of RSTP to GLDP in M2 relative to 
the standard basis X = {x,y}. 
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With each pomt (s, t) G we associate a word 

s t ( —1 — 1\ —t —s 

in F{x,y). Similarly, with a set of points A — {(si, ti), . . . , (s„, i„)} C Z^, 
ordered in an arbitrary way, we associate a word 

n 
1=1 

Observe, that the word Wg^t, as well as wa, belongs to F' = [F, F], so in A/2 they 
define elements from Mj. In particular, the path p^,^ is a closed path in the 
grid T = 1?, which is viewed as the Cayley graph of the abelianization F/F'. 
For A C , (p, (j) e Z^ , and m g Z we put 

• A+{p,q)^{is+p,t + q) I (s,t)eA}, 

• = {(rns,mt) \ {s,t) G A}. 



4 - 

3 r 
2 - 
2 - 
- 
-i- 
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■10 12 3 
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Figure 3: Flows on Z^ defined by words wo,o, W2,i, and u'{(_i.2),(i,-i),(2,3)}- 

The following result is due to Hanan [29] . 

Theorem 4.7. fM Theorem 4J 

Let A — {{xi,yi), . . . , (a;„,?/„)} be a finite subset 0/ Z x Z. There exists 
some Ta for A with the set of Steiner points Q ~ {(oi, &i), ■ • ■ , (a^, bq)} such 
that {fli, C {xi, . . . ,a;„} and {61, . . . ,5^} C {t/i, . . .,y„}. 

Corollary 4.8. Let A be a finite subset of Z x Z. Then 

1. s{A+{b,c)) = s{A) for any (6,c) £ Z^; 

2. s{A) = ms(mA) for any m G N. 

Proof. The first statement is obvious because the parallel shift (x, — > (x, y) + 
(6, c) is an isomorphism of the Cayley graph F (in particular, an isometry). 

To prove the second statement, notice first that s{mA) < ms{A). Indeed, 
stretching Ta by the factor of m along all horizontal and vertical lines gives 
some RST for mA, hence the claim. 
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On the other hand, s{A) < s{mA)/m. To see this, observe that by Theorem 
14.71 there exists R = TmA which hes inside the grid mZ x mZ. Since the 
coordinates of all vertices in R are multiples of k one can shrink R by the 
factor of m, in such a way that the image of R becomes an RST for A. Clearly, 
the size of the image is equal to s(Tm^)/m, hence the result. □ 

Proposition 4.9. Let A be a finite subset of 7?, {b,c) £ A, and n = \A\. Put 
A* = 10n(A- (6,c)). Then Ix{wa') £ [20ns{A),20ns{A) + Aii]. 

Proof. Let m be a geodesic word for wa* relative to the basis X . Since wa* G F' 
the paths p„ and p^,^, are closed paths in T = 1? (viewed as the Cayley graph 
of M2/M2). Hence u and wa* determine the same circulations tt^ = tta^, on 
r. As described in Section [^771 the flow 7r„ is associated with the subgraph r„ 
generated in F by supp(7r„) U {(0,0)}. It follows from the construction of the 
word WA' that the connected components of arc precisely the 1 x 1-squarcs 
in r, whose lower- left corners arc located at the points from A*. Notice that 
(0,0) E A — (6, c) hence (0,0) G A*. Now, if Q is a minimal forest for u (a 
subgraph of F of minimal size that makes the graph r„ U Q connected in F) 
then by Theorem djU] 

\u\=lxiwA')= T^u{e) + 2\E{Q)\=An + 2\E{Q)\. (11) 

Observe, that an optimal RST Ta* for A* also makes the graph r„U(5 connected 
in F, hence \E{Q)\ < s{A*). Therefore, IxiwA') < 20ns{A) + An. 

On the other hand assume that |m| = Ix{wa*) < 20ns{A). Hence, from 
(fTT|) . there exists a minimal forest Q for A* such that 2\E{Q)\ < 20ns{A) — 4n. 
Since every connected component has precisely 4 edges and there are n such 
components, it follows that there is an RST for A* of size strictly less than 
10ns{A) - contradiction with CoroUarv 14.81 This proves the proposition. 

□ 

Corollary 4.10. Let A be a finite subset of and fc G N. In the notation 
above, 

s{A) < k •^=^ IxiwA-) < 2Qnk + 4n. 

In particular, this gives a polynomial reduction of RSTP to BGLP in M2 relative 
to X. 

Proof. Indeed, if s{A) < k then by Proposition [49] /^(wA*) < 20ns{A) + 4n < 
20nfc + 4?!. On the other hand, suppose s{A) > k, say s{A) ^ k + I for some 
positive I e N. Then, again by Proposition l4.9l /y fui4» ) > 20ns{A) ~ 20n{k + 
I) > 20nk + An, as required. □ 

Theorem 4.11. GLDP in a free metabelian group M2 is NP-complete. 

Proof. CoroUarv 14. 101 gives a polynomial reduction of RSTP in 1? to BGLP in 
M2. Therefore BGLP in M2 is NP-liard. Meanwhile, as was mentioned above 
BGLP for M2 is in NP, since WP in M2 is polynomial. □ 
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5 Open Problems 

Denote by Ai the class of all finitely generated metabelian groups. 

Problem 5.1. Describe groups in A4 with GP in P. In particular, the following 
partial questions are of interest here: 

• Are there any groups in M with GP not in P? 

• Do polycyclic groups from A4 have GP in P? 

• Describe wreath products of two f.g. abelian groups have GP in P? 

Notice, since WP is in P for groups from it follows that GLP is at most 
in NP (Turing reducible in P time to BGLP which is in NP). This makes the 
following problem very interesting. 

Problem 5.2. Classify groups in Ai with NP-complete BGLP. In particular, 
the following partial questions are of interest: 

• Do polycyclic groups from Ai have GLP in P? 

• Is it true that a wreath product of finitely generated abelian groups AwrB 
has NP-complete BGLP ii A ^ 1 and the torsion-free rank of B is at least 
2? 

Clearly, GLP is polynomial time reducible to GP. On the other hand, it is 
not clear if there are finitely presented (or finitely generated) groups where GP 
is not polynomial time reducible to GLP. It would be interesting to clarify the 
situation in the class Ai. To this end we post the following 

Problem 5.3. Are there group is A4 where GP is not polinomial time reducible 
to GLP? 
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